First-Order Partial Differential Equations 

Lecture 3 

First-Order Partial Differential Equations 

Text book: Advanced Analytic Methods in Continuum Mathematics, by Hung Cheng (LuBan 
Press, 25 West St. 5th FL, Boston, MA 021 1 1, USA). 

In this lecture and the next two lectures, we'll briefly review partial differential equations 
(PDEs). As PDEs are much more difficult to solve than ODEs, we shall start with the simplest of 
PDEs, those of the first order. 

The good thing about a first-order PDE is this: it can always be "solved" in a closed form. This is 
true whether the PDE is linear or non-linear, and in the former case, whether it is homogeneous or 
inhomogeneous. 

Having made such a strong statement about the first-order PDE, we must add that the solutions 
of the first-order PDE are expressed by the solutions of a system of ODE, which may not be exactly 
solvable. But since ODE is so much simpler to solve than PDE, we consider a PDE "solved" if its 
solution can be expressed as those of an ODE. 

We shall start by making a comparison of the solution of a first-order ODE with those of a 
first-order PDE. 



A Trivial Example 

Let us consider the ODE 

du(x) _ q 
dx 

The answer is, trivially, 

u{x) = c, 

an arbitrary constant. The solution of the ODE is unique if an initial condition is imposed. For 
example, if we require 

u(0) = 1, 

then we have 

In contrast, let us consider the PDE 
where 

The solution of (3.1) is, again trivially, 

u(x,y)=f(y). (3.2) 

Let us describe what this obvious result means in some detail. Eq. (3.2) says that u(x,y) is a function 
of y only, and is independent of x. Therefore, if the value of u(x,y) at an initial point on the y -axis, 
say (0,a), is equal to b, i.e., 

w(0,a) = b, 

then u(x,y) is equal to b on the entire line y = a, i.e., 

u(x,a) = b. 

We call the line y = a a characteristic curve of the PDE (3.1). It is remarkable that u on this 
characteristic curve is independent of the value of u anywhere else. To find the value of u on this 
characteristic curve, we only need the value of u at a point on this characteristic curve, just as with 
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u(x) = 1. 

u x {x,y) = 0, (3.1) 

U x m |«. 

OX 
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the case of the first-order ODE discussed above. Clearly, this initial point does not have to be on the 
y axis. 

If the values of u(x,y) on the y axis between a\ < y < a 2 are given, then the values of w(x,y) are 
known in the strip of the x-y plane with a\ < y < a^. 

If the values of u{x,y) are given on a curve, then the values of u(x,y) are known in a strip 
obtained by sweeping the curve horizontally. 

But if the curve intercepts a horizontal line at two points, then there is a problem. For there can 
be no continuous solution unless the values assigned for u at these two points happen to be equal. 

The situation gets worse if the curve on which the initial data is assigned is a horizontal line. In 
such a case, the initial data is in conflict with the PDE unless all initial values on this horizontal line 
are equal to the same constant. In such a case, the value of the solution at any point not on the 
horizontal line remains unknown. 

While this example is admittedly trivial, the solution of many first-order PDEs can be obtained 
by slightly generalizing what we've just done. 

B Linear Homogeneous PDEs 

Consider the PDE 

a(x,y)u x (x,y) + b(x,y)u y (x,y) = 0, (3.3) 

which is linear and homogeneous. 
Let 

Au(x,y) - u{x + Ax,y + Ay) - u(x,y) 
be the variation of u as the coordinates are varied by (Ax, Ay). By the Taylor expansion, 

Au = u x Ax + u y Ay, 

where the higher order terms, negligible as Ax and Ay become infinitesimal, are ignored. 
Question to the Reader: 
What values should we choose for Ax and Ay in order for Au to vanish? 
Answer 

If u is a solution of (3.3), we have, if a * 0, 

U x = —fiUy. 

Thus 

Au = (— §-Ax + Ay)u y . 
Therefore, Au vanishes if 

Ay = |-Ax, (3.4a) 

or 

—f^T = v^T- (3.4b) 
a(x,y) b(x,y) 

Therefore, u(x,y) is constant on a curve the slope of which is dictated by (3.4b). We may rewrite 
(3.4b) as 

dy _ a(x,y) 
dx b(x,y) ' 

This equation does not depend on u. It is self-contained and can be solved on its own, regardless of 
what u is. We call the solutions of (3.4b) the characteristic curves of the PDE (3.3). 

If a vanishes at some points in the region of interest, there is the difficulty of 1/a blowing up at 
these points. However, if b + at these points, we may express u y by -aujb and get (3.4b). 

If both a and b vanish at a point, then the slope of the charateristic curve may be ill-defined at 
this point. If such is the case, there may be more than one characteristic curves passing through this 
point. 

Problem for the Reader: 
Find the characteristic curves of the PDE 

u t - u x = 0. 

In particular, find the characteristic curve which passes through the point (xo, to). 
Answer 
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The equation for the characteristic curves is 

dt _ dx 
1 1 ' 

the solutions of which are 

a = x + t, (3.5) 

where a is a constant which differs from one characteristic curve to another. In other words, each 
curve is designated by a value of a. Thus the characteristic curves are a family of curves of one 
parameter. 

Conversely, we may easily find the parameter a associated with the characteristic curve passing 
through a point (x, t). For example, the parameter associated with the characteristic curve passing 
through (2.3) is, by (3.5), 

a = 2 + 3 = 5. 

More generally, the parameter a for the characteristic curve passing through (xo, to) is given by 

a = xo + to. 

Problem for the Reader: 
Find the analytic form of the general solution of the PDE 

u t - u x = 0. 

Answer 

Since the values of u are the same on a characteristic curve, its value at (x, t) depends only on the 
value a associated with the characteristic on which the point (x, t) lies. The value a at the point (x, t) 
is equal to (x + t). Thus we have 

u(x,t) = F(x + t) (3.6) 

where F is an arbitrary function. 
Problem for the Reader: 
If the initial values of u(x, t) are given by w(x, 0) = /(x), find u(x, t). 
Answer 

The solution can be expressed in two forms. We may express it in the analytic form by setting 
t = in (3.6). We get 

f(x) = F{x). (3.7) 

Thus 

u{x,t) =f(x + t), 

which is a wave travelling to the left with the uniform velocity unity. 

The solution may also be described graphically by assigning a value of u for each characteristic 
line in the figure on the right. For example, the value of u on the characteristic passing through the 
origin is/(0). 

The analytic solution is simpler to use for the wave equation given above. Indeed, it is simpler to 
use for equations in the form of (3.3). But the graphical solutions are preferable for some of the more 
complicated first-order PDEs not in the form of (3.3). 
Problem for the Reader: 

Find the characteristic curves of the PDE 

u x + (x + y)u y = 0, (3.8) 
and draw some of these characteristic curves in the xv-plane. 
Answer 

The characteristics curves of the PDE (3.8) satisfy 

dx = dy 
1 x + y 

or 

dy 

This is a linear and inhomogeneous ODE of the first order with the complementary solution e x . A 



— 3 — 



First-Order Partial Differential Equations 
particular solution is 

y P {x) = D l _ x = -(1 + D + D 2 + • • »)jc = -x - 1. 

Thus the characteristic curves are given by 

y = ce x - x - 1 , 

or 

e~ x {x +y + 1) = c. 
The characteristic curves for c = 0, 1,-1 are drawn in the figure below. 





/ C= 1, U = 1 


c = D1 , u = e 1 


\ = 0. u = e 
\ i \i 



(3.9) 
(3.10) 



Figure 3.1. 

Problem for the Reader: 
If u{x, y) of (3.8) satisfies the initial condition 

u(0,y) = e-y, (3.11) 
state the values of u(x, y) on the characteristic curves drawn above. 
Answer 

The coordinates of the intercepts of the characteristic curves (3.10) of c = 0, 1, -1 at the y axis 
are y = -1, 0, -2, respectively. According to (3.1 1), the values of u at these three points are e, 1, e 2 
as indicated in the figure above. 

We may obtain u(x, y) in an analytic form. As the value of u depends only on the characteristic 
curve it is on, and as each characteristic curve is defined by the parameter c, u depends only on c; we 
have 

u(x,y) = F(c) = F(e~ x (x + y+ 1)). 
The function F is determined from the initial condition (3.1 1). We have 

e~y = F{y + \), 

or 

F(y) = e l ~y. 

Thus we have 

u{x,y) = exp[l -e~ x (x + y + 1)]. (3.12) 

The solution (3.12) can be viewed geometrically as a surface in the three-dimensional xyu space. 

Non-linear first-order ODE can also be solved in a closed form. These and numerous other 
examples are discussed in the textbook. 
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